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ABSTRACT 

It was found that the equatorial velocity of A-type stars undergoes an acceleration in the 
first third of the main sequence (MS) stage, but the velocity decreases as if the stars were not 
undergoing any redistribution of angular momentum in the external layers in the last stage of the 
MS phase. Our calculations show that the acceleration and the decrease of the equatorial velocity 
can be reproduced by the evolution of the difTerential rotation zero-age MS model with the 
angular momentum transport caused by hydrodynamic instabilities during the MS stage. The 
acceleration results from the fact that the angular momentum stored in the interiors of the stars 
is transported outwards. In the last stage, the core and the radiative envelope are uncoupling, 
and the rotation of the envelope is a quasi-solid rotation; the uncoupling and the expansion 
of the envelope lead to that the decrease of the equatorial velocity approximately follows the 
slope for the change in the equatorial velocity of the model without any redistribution of angular 
momentum. When the fractional age 0.3 < t/tMS ^ 0-5, the equatorial velocity remains almost 
constant for the stars whose central density increases with age in the early stage of the MS phase, 
while the velocity decreases with age for the stars whose central density decreases with age in 
the early stage of the MS phase. 

Subject headings: stars: early-type — stars: rotation — stars: evolution 



1. Introduction 

Helioseismology and asteroseismology have suc- 
cessfully shown us the internal rotation of the 
Sun and red-giant stars. For ex ample, the Sun 
has a flat rotational 



rofile ((Brown et al. 19891: 



Kosovichev et al.Ml997D but red- giant stars have 



a fast rotation core (Aerts et al. 2003: Beck et al 



20121 ). which provided some constraints on the 
processes of angular momentum transport in the 
Sun and red-giant stars. In order to reproduce 
the flat solar rotational profile , the magnetic an- 
gular momentum tran sport (jEggenberger et al 



2005t lYang fc Bil 120061 ) or the gravity-wave an- 



gular momentum transpo rt (jZahn et al.l 11997 



Talon fc Charbonnell 120051 ) should be considered 
in solar models, besides the angular momentum 
transport caused by hydrodynamic instabilities. 
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Most of A-type stars are known to be fast rota- 
tors, wit h a typical rotation al velocity of about 160 



km s ( Rover et al.l 120071 ) for stars with masses 



between about 1.3 and 3.0 Mq. Thus there may 
be a non-negligible effect of rotation on the struc- 
ture and evolution of these stars. However, these 
stars are alway s hotter than the red ed ge of the in- 
stability strip ( Saio fc Gautschv 1998[) . They are 
generally not expected to exhibit solar-like oscilla- 
tions. Therefore the effect of rotation and rotation 
scenarios of these stars are hard to be understood 
as that of the Sun and red-giant stars via seismol- 
ogy- 

Fortunately, the surface rotational veloci- 
ties of thousands of stars with the spectral 
type betw een late B and earl y F have been 
obser v ed (lAbt fc Morrel] 1199,4 IWolff fc Simon 



1997 



Ro ver et al.' '2002 atfhl: IWolff et al I2o1m 

Rover et al. 2004a. b. 2007: lDiaz et al.l201lHZorec fc Rover 
20121) . These stars include pre-main sequence 



(PMS) and main sequence (MS) stars. The ob- 
served surface rotational velocity of a star depends 
on: (1) the initial angular momentum, (2) the 
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rate of angular momentum loss, (3) the angular 
momentum transport in stellar interiors, and (4) 
the changes in the moment of inertia of the star. 
For stars with masses greater than 1.6 Mq, there 
is no firm evidence of activity and any signifi- 
cant angular momentum loss during MS evolu- 
tion (jWolff fc Simonlll997t ). For these MS stars, 
thus, the observed surface rotational velocities are 
mainly dependent on the initial rotational pro- 
file at the zero-age main sequence (ZAMS), the 
angular momentum transport in stellar interiors, 
and the changes in the moment of inertia of stars. 
The changes of the moment of inertia can be un- 
derstood from evolutionary models. Therefore, 
the observed rotational velocities of A-type stars 
could provide important clues about the initial ro- 
tational profile at the ZAMS and the efficiency of 
angular momentum transport in stellar interiors. 

Recently, IZorec fc Roven (|2012l ) studied the 



evolution of surface rotational velocities of more 
than one thousand A-type MS stars. They found 
that the surface velocities of the stars with masses 
between about 1.7 and 3.2 Mq undergo a strong 
acceleration in the first third of the MS evolu- 
tionary phase, which strongly differs from that 
theoretically predicted by two limiting cases of 
internal angular momentum redistributions: (i) 
rigid rotation; (ii) conservation of angular mo- 
mentum in stellar shells, while in the last third of 
the MS the velocities decrease as if the stars were 
not undergoing any redistribution of the angular 
mom entum in the external layers (jZorec fc Rover 

[ioTi. 

The ZAMS models are always assumed to be a 
uniform rotation in the calculations of the evolu- 
tion of rotating stars. If this assumption was cor- 
rect, a special mechanism might be required to ex- 
tract angular momentum from the stellar interiors 
to achieve t he acceleration of the equa torial veloc 



itv. In fact. lWolff fc SimonI (119971 ) and I Wolff et al 



(l2004h lad studied the angular momentum evolu- 
tion of A- and F-type stars from the birth line to 
the MS. They concluded that A- and early F-type 
stars do not loss angular momentum and the an- 
gular momentum could be conserved in the shells 
of the stars as the stars evolve from the end of 
fully convective phase to the ZAMS. Thus, for 
these stars the rotation at the ZAMS should be dif- 
ferential (hereafter this differential rotation 
ZAMS model obtained from the evolution 



with the conservation of angular momen- 
tum in shells was referred to Wolff ZAMS 
model), and the core could serve as a reservoir 
of angular momentum which can be transported 
from the fast rotation core to the slow rotation 
envelope to produce the acceleration of the equa- 
torial velocity. Moreover, the too high or too low 
efficiency of angular momentum transport in stel- 
lar interiors could not reproduce the acceleration 
of the surface veloci ty in the first third of the MS 
evolutionary phase dZ orec fc Roved [20i3). Thus 
the results of IZorec fc Rover (2012) p rovide us an 



oppor tunity to test the conclusions of I Wolff et al 



(|2004[ ) and the angular momentum transport in 
the interiors of A-type stars. 

In this work, we focus mainly on the evolution 
of the equatorial velocity of A-type stars. The pa- 
per is organized as follow. We simply show our 
stellar models in section 2. The results are repre- 
sented in section 3. In section 4, we discuss and 
summarize the results. 

2. Stellar models and calculation results 
2.1. Stellar models 



We 



used the Yal e Ro tation Evolution Code 
19891: 



( Pinsonneault et al 



Yang fc Bil 12007 ) to 



compute the evolutions of rotating models with 
different masses and metallicities. Hydrody- 
namical instabilities considered in this code are 
meridional circulation, the Goldreich-Schubert- 
Fricke instability, and the secular shear ins t abilit y 
(jEndal fc Sofial Il978l : IPinsonneault et all [l989 ). 
The OPAL EOS tables (iRogers fc Navfonov|2no; ), 
OPAL opacity tables (jlglesias fc Roger si Il99(j) . 
and the opacity tables for low tempe rature pro- 
vided bv [Alexander fc Ferguson (1994) were used. 
Energy transfer by convection is treated according 
to the standard mixing length theory. The value 
of 1.72 for the mixing-length parameter (a) was 
calibrated against the Sun. The initial metal- 
licity Z was fixed at 0.02 and 0.008, and 
the hydrogen abundance was determined by 
X = 0.767 — 3Z. For a given mass, the initial an- 
gu lar momentum i s estimated by using the formula 
of iKawalen (jl987l ). Angular momentum loss due 
to magnetic braking and mass loss could be negli- 
gible in A-type stars ( MacGregor fc Charbqnneau 



1991 IWolff fc Simonlll997t IZorec fc Roveill2012l ) 

Thus we assumed that the total angular momen- 
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turn is conserved in our models. 

We calculated the following four evolutionary 
cases: (1) CI, models were evolved from the PMS 
to the terminal age main sequence (TAMS) with- 
out any exchange of angular momentum in shells; 

(2) C2, models were computed as a rigid rotator; 

(3) C3; and (4) C4. In the last two cases, mod- 
els were evolved from the ZAMS to the TAMS 
with angular momentum transport caused by the 
rotationally induced instabilities, and the ZAMS 
model is a uniform rotator for the C3 but is a 
Wolff ZAMS model for the C4. Moreover, we 
assumed that the rotation of the surface layers (the 
fractional mass 5M/M ~ 10"'*) is uniform. 

2.2. Calculation results 

Figure [T] shows the evolution of the equatorial 
velocity of models with M = 1.7 M© and Z = 0.02 
in the MS life time span (tMs)- The notation 
VzAMS is the equatorial velocity at the ZAMS. 
For the cases CI, C2, and C3, from the ZAMS to 
t/tzAMS ^ 0.95, the equatorial velocities decrease. 
In case CI, the equatorial velocity descends about 
45%, which is consistent with the expectation of 
the law of conservation of angular momentum in 
shells {Ve = VoRo/Re, where Re is the radius of 
stars and Vq and Rq refer to initial values) . In the 
case C2, the equatorial velocity decreases around 
20%. These two limiting cases did not reproduce 
the acceleration of the surface velocity in the early 
stage of MS. The evolution of the equatorial veloc- 
ity of case C3 is similar to that of the case C2 ex- 
cept for in the last stage of the MS. This is because 
the rotation of the radiative region of the MS mod- 
els of the C3 is a quasi-solid body rotation and the 
core and the envelope are uncoupling in the last 
stage. In the case C4, the equatorial velocity un- 
dergoes an acceleration in the early stage of the 
MS phase. From t/tus - 0.017 to t/tus - 0.3, 
the equatorial velocity increases about 20%. But 
when t/tMS increases from about 0.3 to 0.5, the 
velocity remains almost constant. After t/tus > 
0.5, the velocity decreases as evolution proceeds 
until the evolution approaches the end of the MS 
stage. When the fractional age t/tus > 0.7, the 
decrease of the equatorial velocity of the case C4 is 
similar to that of t he case CI, which is c onsistent 
with the finding of IZorec fc Roveil (|2012[ ). 

Figure [2] shows the radial distributions of in- 
ternal angular velocity of the model with M = 



1.7 Mq at different stages of the evolution of case 
C4. Because there is no redistribution of angu- 
lar momentum in shells during the PMS stage 
and the core of the model contracts but the en- 
velope expands, the core rotates faster than the 
envelope when the model evolves to the ZAMS. 
From t/tus — 0.017 to t/tus — 0.30, the angu- 
lar momentum stored in the fast rotation core is 
transported to the slow rotation envelope, which 
almost fully impedes the decrease in the angular 
velocity of the envelope, even leads to an increase 
in the angular velocity. In addition, due to the in- 
crease in radius, the equatorial velocity increases 
obviously. When the model evolves from t/tus ~ 
0.3 to t/tMS ~ 0.5, the core and the envelope 
are strongly coupled by hydrodynamic instabili- 
ties. Although the core still contracts, its angu- 
lar velocity decreases with the decrease in the an- 
gular velocity of the envelope. The angular mo- 
mentum stored in the stellar interiors is contin- 
ually transported outwards, which is insufficient 
to remain the acceleration of the equatorial veloc- 
ity but can allow the velocity to remain almost 
constant. When t/tus > 0.6, due to the facts 
that the core contracts and the envelope expands 
rapidly, the hydrodynamical instabilities are insuf- 
ficient to couple the core to the envelope. Thus the 
core and the envelope are uncoupling. However, 
the instabilities are sufficient to keep the enve- 
lope rotating as a quasi-solid body. In our models, 
the differential rotation results from the contrac- 
tion/expansion of stars. Due to the uncoupling of 
the core and the envelope, the decrease of the an- 
gular velocity of the envelope depends almost only 
on the expansion of the envelope. Moreover, Fig. |3] 
shows that most of the radiative region spin down 
at an approximately equal rate from t/tMS ~ 0.7 
to t/tus ~ 0.9 in the evolution of case CI, which 
is similar to the manner of quasi-solid rotation. 
Thus when the fractional age t/tus > 0.7, the de- 
crease in the equatorial velocities of cases CI and 
C4 follows an approximate slope. 

We also calculated the evolutions of the 
equatorial velocity of models with M = 2.0, 
2.5, and 3.0 M0. The comparisons of the 
evolutions of the equatorial ve locities of our 
mode ls with the results of iZorec &: Rover 
(|2012( ) are shown in Fig. |4l The evolutions of 
cases CI, C2, and C3 of these models are similar to 
those of model with M = 1.7 Mq, i.e. the equa- 
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torial velocity decreases with age. But in the 
early stage (t/tMS ^ 0.3) of the MS of the case C4, 
the equatorial velocity undergoes an acceleration. 
However, the magnitude of the increase of the 
ratio Ve/VzAMS of our models is obvio usly lower 
than that obtained bv lZorec &: Rove r (2012). 
For example, at t/tus = 0.35, the value of 
the observed Vc/Vzams is 1.265 ± 0.079 for 
the star with AI = 2.5 Mq, but the value for 
our models is about 1.11. The discrepancy 
is about 2a. When the fractional age t/tus in- 
creases from around 0.3 to 0.5, the theoretical 
velocities decrease slightly but the observed 
ones decrease rapidly for stars with AI — 2.5 
and 3.0 Mq. When the age t/tus > 0.5, the 
value of the theoretical Ve/VzAMS is lower 
for the star with AI — 2.0 Mq but is higher 
for the star with AI = 3.0 Mq than the ob- 
served one. However, for the star with 
AI = 2.5 Mq, the evolution of the theoret- 
ical Vo/VzAMS is good agreement with the 
observed one. After t/tMS > 0.7, the decrease 
of the velocities of cases CI and C4 and the ob- 
served one follow an approximate slope. 

Furthermore, the Fig. [4] shows that 
the evolution of the equatorial velocities 
is slightly affected by the metallicity. This 
can be due to the fact that the radius of 
low-Z model is smaller than that of high-Z 
model. 

3. Discussion and Conclusions 

The increase of the fractional velocity Vc/Vzams 
is mainly dependent on the angular momentum 
that is stored in the interiors of the star and can 
be transported outwards to accelerate the surface 
layers. The angular momentum depends on the 
contraction of the stellar core and the expansion 
of the envelope of the star in our models. The 
more massive the star, the less the contraction of 
its core when the star evolves from the PMS to 
the ZAMS. In addition, in the early stage of the 
MS, the central density increases as the evolution 
proceeds for the stars with AI < 2.0 Mq, but the 
central density decreases (see the Fig. [5]) for the 
stars with AI > 2.1 Mq, i.e. the angular veloc- 
ity of the core decreases for stars with AI > 2.1 
Mq. As a consequence, the more massive the star, 
the smaller the fractional angular momentum that 



can be transported outwards. Moreover, the more 
massive the star, the larger the expansion of the 
radius during the MS stage, thus the more difficult 
to remain the angular velocity of the surface layer. 
Therefore the increase of the fractional velocity de- 
creases with increasing the mass, and the available 
angular momentum stored in the stellar interiors 
is almost drained when t/tus ~ 0.3 for the stars 
with M > 2.1 Mq. Hence from t/tMS ~ 0.3 to 
t/tus ~ 0.5, the equatorial velocity decreases for 
stars with M > 2.1 Mq. 

The instabilities caused that the angu- 
lar momentum was transported outwards 
in the early stage of the MS, which leads 
to the increase in the Vo/Vzams- Although 
the magnitude of the increase in the fractional 
vel ocity of o ur models is less than that obtained 
by IZorec fc R over (2012). who given the magni- 
tude as high as 30%, the velocity increases as the 
evolution proceeds during the first third of the 
M S phase, which is co nsistent with the finding 
of IZorec fc Roveil (|2012l) . This might offer a sup- 
port for the Wolff ZAMS model and imply 
that the core and the envelope is uncou- 
pling as A-type stars evolve from the end 
of fully convective phase to the ZAMS but 
is coupling in the early stage of the MS. The 
theoretical magnitude is lower than the ob- 
served one. This might be caused by that 
the efficiency of angular momentum trans- 
port wasn't calibrated. 

In the last stage of the MS (0.7 < t/tus < 0.95), 
the core and the radiative envelope are decoupling 
in the evolution of cases CI, C3, and C4. The 
contraction of the core almost not affect the rota- 
tion of the envelope. For the case CI, although the 
rotation of the whole radiative envelope is differen- 
tial, the envelope spins down at an approximately 
equal rate in our models, which is similar to the 
manner of the quasi-solid body rotation. For the 
case C4, the rotation of the radiative envelope is 
a quasi-solid rotation. Thus the decrease of the 
fractional velocity of the cases CI and C4 follows 
an approximate slope when t/tus ^ 0.7. As a con- 
sequence, although the velocities of A-type stars 
decrease as if the stars were not undergoing any 
redistribution of the angular momentum in their 
envelope during the final stages of the MS phase, 
the rotation of the envelope of A-type stars might 
be a quasi-solid rotation with the decoupling of 
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the core and the envelope. Thus, the rotation 
of A-type stars might be closer to the quasi- 
solid body rotation during t/tus ~ 0.3 - 0.5 
than in the early or last stage of the MS. 
Compared with lEk strom et al. ('2008') 



models used by IZorec &: Rover ( 2012 ), our 
ZAMS models for C4 are the Wolff ZAMS 
model but the rotation of the ZAMS mod- 
els of Ekstrom e t al. (2008) and our C3 is a 
solid body rotation. The equatorial veloc- 
ity of our C3 decreases wit h age, which is 
similar to th at obtained by IZorec &: Rover 
(|2012h from lEkstrom et alJ (|2008h models. 
The cases C3 and C4 have the same in- 
stabilities, this may imply that the Wolff 
ZAMS model is required to achieve the ac- 
celeration. Moreover, we neglected the an- 
gular momentum loss caused by mass loss. 
Finally, there are some differences in the 
treatments of instabilities included models 
([Ping o nneault et al. 1989; Maeder &: Zah^ 
1998h . which can lead to the differences in 
the efficiency of angular momentum trans- 
port and the mixing of elements. Besides 
the hydrostatic and Von Zeipel effects, the 
mixing caused by rotationally induced in- 
stabilities is one of the key factors that af- 
fects the distribution of density and the 
instability of convection by changing the 
distribution of elements. In our models, 
the effects of rotation lead to a decrease 
in the convective core for stars with M < 
2.0 M© but an increase in the convective 
core for stars with M > 2.1 Mq (jYang et al. 
2012 h. As a consequence, near the end of 



the MS, the radius of rotating models is 
smaller than that of non-rotating ones for 
the stars with M < 2.0 Mq but larger than 
that of non-rotating models for the stars 
with M > 2.1 Mq. Due to the fact that the 
decrease of the ratio Vc/Vzams is relevant to 
the increase in the radius, thus the instabil- 
ities make the stars with M < 2.0 Mq easier 
remain their Vo/Vzams at a high value than 
the stars with M > 2.1 Mq. 

In this work, we calculated the evolution of 
the equatorial velocity of A-type stars with dif- 
ferent conditions. The evolution of the case C4 
reproduces the most characteristics of the equa- 
torial velocity of A-type stars that were found 



by IZorec k Roveil (|2012[ ). From the ZAMS to 
t/tMS ~ 0.3, the equatorial velocity undergoes an 
acceleration. But in the last stage of the MS stage 
(t/tMS ^ 0.7), the core and the envelope of A- 
type stars are uncoupling; the envelope rotates as 
a quasi-solid body, but the decrease of the equato- 
rial velocity approximately follows the slope of the 
velocity of the model without any redistribution 
of the angular momentum in this stage. However, 
the magnitude of the increase of the fractional ve- 
locity Vn/VzAMR of our models i s less than that 
obtained by Zorec fc Roveij (l2012h . Moreover, our 
models with M = 2.5 and 3.0 Mq did not repro- 
duce the rapid decrease in the equatorial velocity 
during the middle stage of the MS phase. This 
might imply that some mechanisms for the 
transport or loss of angular momentum are 
missing in our models. 
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Fig. 1. — Evolution of the equatorial velocity 
in the MS stage. The long-dashed (red) line 
shows the evolution of CI. The dash-triple-dotted 
(green) line indicates the result of C2. The solid 
(blue) line is given for C3. The dotted (cyan) line 
corresponds to the evolution of C4. 
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Fig. 3. — Same as Fig. [5] but for the evolutions of 
case CI. 
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Fig. 2. — The radial distributions of internal 
angular velocity at different evolutionary stages. 
The solid (red) line shows the distribution at 
Xc = 0.703. The long-dashed (green) line indi- 
cates the distribution at X^. — 0.535. The dash- 
triple-dotted (blue) line shows the distribution at 
X(, — 0.451, and the dotted (cyan) line refers to 
the distribution at Xr = 0.254. 




Fig. 4. — Comparisons of the evolutions of the- 
oretically equatorial velocities with the observed 
data. T he dash-dotted li nes (g reen) show the re- 
sults of IZorec fc Rover ( 2012 ). The error bars 
(green) represent Icr errors. 
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Fig. 5. — The central density (pc) as a function of 
the mass fraction of central hydrogen. The dashed 
(red) lines show the results of models without ro- 
tation. The dotted (cyan) lines indicate the results 
of models with rotation. 
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